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Among various strong-curvature extensions of general relativity, Einstein-dilaton-Gauss-Bonnet 
gravity stands out as the only nontrivial theory containing quadratic curvature corrections while 
being free from the Ostrogradsky instability to any order in the coupling parameter. We derive an 
approximate stationary and axisymmetric black hole solution of this gravitational theory in closed 
form, which is of fifth order in the black hole spin and of seventh order in the coupling parameter of 
the theory. This extends previous work that obtained the corrections to the metric only to second 
order in the spin and at the leading order in the coupling parameter, and allows us to consider values 
of the coupling parameter close to the maximum permitted by theoretical constraints. We compute 
some quantities which characterize this solution, such as the dilaton charge, the moment of inertia 
and the quadrupole moment, and its geodesic structure, including the innermost stable circular 
orbit and the epicyclic frequencies for massive particles. The latter provides a valuable tool to test 
general relativity against strong-curvature corrections through observations of the electromagnetic 
spectrum of accreting black holes. 


I. INTRODUCTION 

Observational and experimental tests of general 
relativity (GR) [T] have mostly probed the weak- 
field/slow-motion regimes of the theory, while a 
number of strong-held, relativistic GR predictions 
still remain elusive and difficult to verify [1 E]. 
Furthermore, a series of long-lasting problems in 
Einstein’s theory - such as the accelerated ex¬ 
pansion of the Universe, dark matter, the nature 
of curvature singularities, and the quest for an 
ultraviolet completion of GR - have motivated 
strong efforts to develop extended theories of grav¬ 
ity which would modify GR in its most extreme 
regimes while conforming with current weak-held 
observations [3]. 

Black holes (BHs) are genuine strong-held pre¬ 
dictions of GR and have no analog in Newto¬ 
nian theory. Thus, they are natural candidates 
to test gravity in the strong-held regime. Future 
networks of electromagnetic detectors [MS] and 
ground-based gravitational-wave detectors ladni 
will allow us to measure some crucial properties of 
BHs, such as their shadows, the location of the 
event horizon, and of the innermost stable cir¬ 
cular orbit (ISCO). This information will be in¬ 
strumental to test the Kerr hypothesis, accord¬ 
ing to which all stationary astrophysical BHs are 
uniquely described by the Kerr family and are, 
thus, characterized by only two parameters: their 
mass and angular momentum (see, e.g.. Ref. m 
and references therein). 

In recent years, several modihed theories of grav¬ 
ity have been proposed. They can be divided in 


various categories, each one lifting some of the 
fundamental principles (Lorentz invariance, weak 
and strong equivalence principles, massless spin-2 
mediators, etc.) upon which Einstein’s theory is 
uniquely built [3]. From this and other classifica¬ 
tions, it emerges that one of the simplest and best 
motivated ways to modify GR consists of including 
a fundamental scalar field which is nonminimally 
coupled to the metric tensor. In order to modify 
the strong-curvature regime, it is natural to couple 
this scalar field in the gravitational action to terms 
quadratic in the curvature tensor. Such couplings 
can also be interpreted as the first terms in the ex¬ 
pansion in all possible curvature invariants, as sug¬ 
gested by low-energy effective string theories [Hj. 
Generally, a quadratic curvature term in the ac¬ 
tion leads to field equations of third (or higher) or¬ 
der, which are subject to Ostrogradsky’s instabil¬ 
ity m- Therefore, these theories should be con¬ 
sidered as effective, i.e., truncations of a theory 
with further terms in the action, which are neg¬ 
lected in the perturbative regime. 

It should also be mentioned that quadratic 
curvature terms are crucial, not only to modify the 
strong-curvature regime of GR, but also to affect 
the behavior of stationary BHs. Indeed, standard 
scalar-tensor theories (in which one or more scalar 
fields are included in the gravitational sector of the 
action) satisfy the so-called no-hair theorems; i.e., 
stationary, vacuum BHs are the same as in GR 
[T3HIS] (but see UMISj for possible violations of 
these theorems). When quadratic curvature terms 
are included in the action instead, stationary BH 
solutions are different. 


2 


We shall consider a member of this family of 
modified-gravity theories, Einstein-dilaton-Gauss- 
Bonnet (EDGE) theory, in which a scalar field 
(the dilaton) is coupled to the Gauss-Bonnet in¬ 
variant [ 121120 ] 

^GB = + R^ (1) 

in the action. EDGE gravity is one of the best 
motivated alternatives to GR. Indeed, it is the 
only theory of gravity with quadratic curvature 
terms in the action, whose field equations are of 
second differential order for any coupling, and not 
just in the weak-coupling limit which is assumed 
in the effective-field-theory approach |4]. As a 
consequence, EDGE gravity is ghost-free; i.e., it 
avoids the Ostrogradsky instability m- Eurther- 
more, as mentioned above, the higher-curvature 
coupling - which modifies the strong-curvature re¬ 
gime of gravity - violates the hypothesis of the EH 
no-hair theorems so that EH solutions in EDGE 
gravity are different from those predicted by GR 
and provide the ideal arena for genuine strong-held 
tests of the Kerr hypothesis. Finally, the EDGE 
term naturally arises in low-energy effective string 
theories [21] . 

In this work, we construct an analytical, per¬ 
turbative solution of EDGE theory, which de¬ 
scribes a slowly rotating EH endowed with a scalar 
held. To this aim, we extend the formalism de¬ 
veloped in [221 123] up to hfth order in the EH 
(dimensionless) spin parameter y = where 

J and M are the angular momentum and the 
Arnowitt-Deser-Misner mass of the solution, re¬ 
spectively. 

Analytical EH solutions of EDGE theory in 
the small-coupling limit have been investigated 
in [MIES], where stationary, spherically symmet¬ 
ric conhgurations where found. Approximate, sta¬ 
tionary, and axisymmetric solutions to linear and 
quadratic order in the EH spin were obtained in 
[2S] and [27], respectively. Both of these works 
considered a weak-held expansion of the coupling 
between the scalar held and the Gauss-Bonnet in¬ 
variant in terms of a dimensionful coupling 

constant a. Exact numerical solutions were con¬ 
structed to zeroth [20] and hrst order [2H] in the 
spin and also for arbitrary values of the angular 
momentum [291 Although exact in a, such 
solutions are of limited practical use (for instance, 
for Monte Garlo data analysis) because they re¬ 
quire a numerical integration for each set of para¬ 
meters. On the other hand, numerical solutions are 
necessary in regimes where the slow-spin expansion 
does not converge and are, therefore, complement¬ 
ary to our analysis. 

Our results extend the study carried out so 
far. In particular, we go beyond the analysis of 
Ref. [27], where a EH solution was obtained to 


second order both in the spin and in the coup¬ 
ling parameter. Indeed, we compute the metric 
tensor and the scalar held up to O where 

C = a/M^, and a is the EDGE coupling constant. 
We use this expansion to derive the main features 
of the solution, such as the geometry of the event 
horizon and of the ergoregion. Furthermore, we 
study the geodesic structure of this solution by 
computing the IS GO and the epicyclic frequencies 
(see, e.g.. Refs. [3TH53] 1 consistent with our ap¬ 
proximation scheme. We compare these quantit¬ 
ies with those obtained in |h4| . where a numerical 
solution was derived, which is exact in the coupling 
parameter (i.e., with no perturbative expansion in 
C) and approximate to linear order in the BH spin. 
We find relative errors at most of the order of 1% 
for the maximum value of ^ allowed by theoret¬ 
ical constraints for the existence of BH solutions 
C ^ 0.691 [20] and much smaller for less extreme 
couplings. 

The results of this paper can be useful to devise 
tests of GR in the strong-held regime through as- 
trophysical observations of BHs. For instance, we 
have shown |34| that observations of quasi-periodic 
oscillations of accreting BHs, with the sensitivity 
of recently proposed large-area X-ray space tele¬ 
scopes (e.g., 017]), allow us to set constraints on 
the parameter space of EDGE theory, thus, prob¬ 
ing the strong-held regime of gravity (see, also. 
Ref. [22] for a recent study). However, since BH 
solutions in EDGE theory (for hnite a) were only 
known at hrst ordeiQ in the spin parameter y, 
in [34] we only considered BHs with very slow- 
rotation rate, for which the deviations from GR 
are expected to be small. 


This paper is organized as follows. In Sec. [TT] 
we derive our solution of the EDGE held equa¬ 
tions, describing rotating BHs up to C>(C^,y®). 
In Sec. jlllj we study this solution, computing its 
geometrical properties, the location of the ISGO, 
and the azimuthal and epicyclic frequencies. We 
also estimate the accuracy of our approximation 
in the determination of these quantities and how 
our results improve on the existing literature. In 
particular, we discuss how the spin correction to 
the azimuthal and epicyclic frequencies can affect 
possible tests of GR based on observations of ac¬ 
creting BHs, such as those discussed in [23] ■ Fi¬ 
nally, in Sec. IV we draw our conclusions. 


^ As mentioned above, a solution for finite spin and coup¬ 
ling is only known in numerical form 12911301 , and it is 
impractical for extensive studies of geodesic properties. 
However, numerical solutions are necessary to explore the 
high-spin regime, especially because EDGB BHs can vi¬ 
olate the Kerr bound and can have x > 1 m- 
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II. SPINNING BLACK HOLES IN 
EINSTEIN-DILATON-GAUSS-BONNET 
THEORY 

In this section we derive the spacetime metric 
and scalar field describing rotating BHs in EDGE 
theory, up to 0 {(^,x^)- 


tensor where 

- a/C^, , ( 6 ) 

is the effective stress-energy tensor for the dilaton. 

B. Static BH solutions 


A. EDGB gravity 


Einstein-dilaton-Gauss-Bonnet theory is defined 
by the following action [121 EQ] : 


1 OP* 

( 2 ) 

where g < 0 is the metric determinant, d* is a 
scalar field coupled to the Gauss-Bonnet invari¬ 
ant 0 and a > 0 is the coupling constant |5D]. 
Since we are interested in BH solutions, in the ac¬ 
tion above we have neglected matter fields. We 
use geometric units G = c = I: with this choice, 
the scalar field $ is dimensionless and a has the 
dimensions of a length squared. 

The field equations of EDGB gravity are found 
by varying the action ([^ with respect to and 



G^^, = - alC^, , (3) 

5 = = 0, (4) 

where G^i, = — \gfi.uR is the Einstein tensor, 

\ i9i^p9<y\ + 9pX9up) 

xV^ , (5) 

and is Levi-Givita tensor, with _ 

— (— 5 )“^/^. Note that - by virtue of the GB 
combination entering the action 0 - the equa¬ 
tions are of second differential order, and, there¬ 
fore, this theory is free from the Ostrogradsky in¬ 
stability [13]. Indeed, EDGB gravity is a particular 
case [36] of Horndeski gravity - the most general 
scalar-tensor theory with second-order field equa¬ 
tions m- This special subcase is the only one 
known to date in which regular, stationary, asymp¬ 
totically flat, hairy BH solutions other than GR 
ones are found [5B]. Furthermore, EDGB grav¬ 
ity can be obtained from the low-energy expan¬ 
sion of the bosonic sector of heterotic string the¬ 
ory [H] |3T| , in such case the coupling a is related 
to the string tension. 

In order to simplify our notation, in the next 
sections we shall introduce the modified Einstein 


Since the EDGB coupling constant has the di¬ 
mensions of the inverse of the curvature tensor, 
it is natural to expect that in this theory the 
strongest deviations from GR will come from phys¬ 
ical systems involving high curvature, such as BHs, 
neutron stars and the early Universe. We focus 
here on BH solutions and, in particular, on ro¬ 
tating BH geometries that are obtained through 
a slow-rotation expansion around a static back¬ 
ground solution. 

The exact BH background solution (first derived 
in [ 20 ]) is described by the static, spherically sym¬ 
metric line element 

ds^ = + r^dn^ , (7) 


and by a spherically symmetric scalar field, $ = 
^(r). The field equations (|^ and Q supplied by 
the metric ansatz 0 reduce to a set of differential 
equations for the scalar field and for the functions 
r and A. Indeed, Eq. Q yields 



+ (l-e-^) 


r' 

r" + _(r' 



, ( 8 ) 


while the t-t, r-r and 9-9 components of G^i, = 0 
reduce to 


1 + _^'(i _ 3e-A) 


A' = 




r I — 


ae^ 


4 r 

I il2\ 


r 


1 + ^^'(1 _ 3 e -^) 


r' 1 


r' = 


-k'2 


r e^-l 


( 10 ) 


2 r 


- + - (r'-A') = -V 


(j)'^ ae' 


<f>-A 


2 

r'x/ 


^T" + r'((^" + <^'2) + _ 3A') 


( 11 ) 


Note that Eqs. 0 ~( 1 TT]) are not all independent 
and that the r-r component can be solved analyt¬ 
ically, yielding 


e^ = 


-f3 -b a//32 - 47 


( 12 ) 


2 
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where 


C. Spinning BH solutions 


/3 = 


4 


-1-r' 



7 = 2r>'e^ 


(13) 

The remaining two independent equations can be 
written as 



where the radial functions d, di and d 2 are given 
in Appendix A of [5D] . The Arnowitt-Deser-Misner 
mass M and the dilatonic charge T) can be read off 
the asymptotic behavior of the metric and of the 
dilaton field, 


To describe slowly rotating BH solutions we ex¬ 
tend the approach developed by Hartle [221 US] j in 
which spin corrections to the static solutions are 
introduced within a perturbative framework. The 
procedure described in this section is generic and 
can be applied also to other theories and different 
spinning solutions. 

Let us start with the most general solution for a 
stationary, axially symmetric spacetime]^ which is 
given by 

ds^ = -H^df + Q^dr^ + r^K^[d9^ + 

+ sin^ 0{d(p — Ldt)^] , (19) 


Qtt ~1 + ‘i.M/r -f ... (15) 

(j) V jr -f ... (16) 


It turns out that for each value of M, there is only 
one solution describing a static BH. In other words, 
the scalar field is a “secondary hair”: the dilatonic 
charge T) is not an independent parameter but is 
determined in terms of the BH mass M. 

The field equations are invariant under the res¬ 
caling (j) ^ (j) + ^ and r —>■ (or, equivalently, 

M —>■ and T> where 0 is a con¬ 

stant. We fix this freedom by requiring that 0 
at spatial infinity; this means that at infinity, the 
Gauss-Bonnet invariant appears in the action ([^ 
multiplied by the constant a/4. 

As noted in [20] static BH solutions in EDGB 
gravity exist only if 




r\, 

a-\/6 


(17) 


where 0h is the value of the scalar field computed 
at the horizon rh. As shown in [28j . by requiring 
that (/ —)■ 0 at spatial infinity, Eq. ( [Tt] ) can be 
recast in the form 


where H,Q,K, and L are functions of {r,0). The 
ansatz (191 can be expanded perturbatively in the 
spin around the static solution 


ds^ =—e'"[l + 2h{r, 9)]dt^ + e ^[1 -I- 2m(r, 6)]dr^ 

-|-r^[l -I- 2fc(r, 9)][d9^ -\- sin^ 9{dip — uj{r, 9)dt )'^], 

( 20 ) 


where the functions w, h, m, and k can be expan¬ 
ded in a complete basis of orthogonal functions ac¬ 
cording to their symmetry properties as 



^x-<? 

n 

xM'^\r)Sii9), 


UJ = 

E 

E 

(21) 


n=l,3,5,. 





Nx-P 

n 

X^h\^\r)Pi {cos 9) , 


h = 

E 

E 

(22) 


n=2,4,... 

/^0,2,4,... 




Nx-P 

n 

x'^m^P\r)Pi {cos 9) , 


m = 

E 

E 

(23) 


n=2,4,... 

Z=0,2,4,... 




Nx-P 

n 

x"’k\'^\r)Pi{cos9) , 


k = 

E 

E 

(24) 


n=2,4,... 

Z=0,2,4,... 




0<^<0.691. 


(18) 


Thus, smaller BHs would correspond to a more 
stringent bound on a. 

Presently, the tightest observational bound on 
the EDGB coupling parameter (obtained by the 
orbital decay of X-ray binaries) is a < 47M|| [55] . 
As discussed in [HI], this upper bound is weaker 
than the theoretical constraint 0 for BHs with 
M < 8.2Mq. Eor such BHs, the entire range (181 
is phenomenologically allowed. 


Solutions of Eqs. (14) have been solved numer¬ 
ically in Ref. [20] , whi e an analytical static BH 
solution has been derived to second order in a/M^ 
in Refs. 1551 [551. 


where Pi are the Legendre polynomials. Si = 

Po = Si = 1), and p 
(respectively q) is zero when the order of the 
spin expansion is even (respectively odd), whereas 
p = 1 (resp. g = 1) otherwise. The radial 

functions the order 


©(x"). Note that, since the metric (20) is invari¬ 
ant under the rescaling r —>■ /(r), the functions 
k^\r) can be set to zero without loss of general¬ 
ity [21[231- 


^ We also assume equatorial symmetry and invariance un¬ 
der (t —t, ip —V?)- 
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Because the dilaton field transforms as a scalar 
under rotation, we expand it as 

N^-p n 

$(r) = (/)(r)+ ^ ^ x”4”^(?')^i(cos6») , 

n=2,4,... Z=0.2,4,... 

(25) 

where 4 > is the background static solution and the 
radial functions are of the order 0(x”). 


1. 0{x) corrections 


Rotating BH solutions in EDGB gravity have 
been investigated to linear order in the spin an¬ 
gular momentum in Refs. [5S1[55]. At this order, 
the metric (20) reduces to the static case with a 
nonvanishing gravitomagnetic term described by 
a)(r, 0) = = uj{r) [see Eq. ([^]: 

ds^ = + r'^dVL^ 

—2r^a;(r) sin^ Odtdip . (26) 


Erom = 0, it is easy to show that uj satisfies 
the second-order equation ^5] : 

oj”+ [2r^e^ — 2ar(/)'e‘^] 

(-ae"^ [ 2 (i)''r + (()'(6 -h 20'r - E'r - 3AV)] 

-e^r[-8-hr(r'-hA')])a;' = 0, (27) 

where the coefficient of ui' depends on the nonspin¬ 
ning solution. The BH angular momentum can be 
read off the asymptotic behavior of the gyromag- 
netic term. 


tLi(r) —>■ 


2J 


(28) 


at large distance. 


2. 0{x") corrections: n > 2 and even 


Replacing the metric ansatz (20) into the field 
equations and using the decomposition in Legendre 
polynomials, a set of ordinary differential equa¬ 
tions can be obtained, at each order in the spin 
expansion. The equations are inhomogeneous with 
source terms given by the lower-order functions. 

At each given order n > 2 with even n, the equa¬ 
tions are found from Ei = Gtt = 0 , E 2 = Grr = 0, 
E 3 = Gee + {sin9)~'^Gtf,tf, = 0, and if 4 = 5 for 
the scalar equation , each contracted with a Le¬ 
gendre polynomial, 


d9sin6Pi{cos9)Ei{r,6) = 0 , (29) 

where i = 1,2,3,4 and I = 0, 2,4,..., n. Because 
of the symmetry properties of the field equations 



and of the background, this procedure gives a set of 
purely radial, inhomogeneous, ordinary differential 


equations for m 
0,2,4,... ,n (we recall that = 0). 


and 


with I = 


3. G(x") corrections: n > 3 and odd 

Similarly, at a given order n > 3 (with odd n) 
in the spin expansion, a set of radial equations for 
the gravitomagnetic terms can be obtained by con¬ 
tracting Gtip = 0 with the (axisymmetric) vector 
spherical harmonics, namely 

d0sin0^^^^^°^Gt^ = O (30) 
d cos 9 

with I = 1,3, 5,..., n. Again, this procedure yields 
a set of purely radial, inhomogeneous, ordinary dif¬ 
ferential equations for with / = 1, 3, 5,..., n. 



D. Small-coupling approximation 


The set of equations presented above provides 
a full description of the BH solution at any per¬ 
turbative order in the spin, but generic (i.e., non- 
perturbative) in the EDGB coupling. However, 
such equations are cumbersome and it is imprac¬ 
tical to solve them numerically. More importantly, 
the theoretical constraint (18) shows that the di¬ 
mensionless coupling parameter has to be smaller 
than unity. This motivates a small-coupling ap¬ 
proximation pSll28j . in which the field equations 
are solved perturbatively in a/M"^ ^ 1 to some 
desired order. Actually, because we are interested 
in the regime a/M^ < 1 (the maximum valu^of 
this parameter is 0.691), we shall compute terms 
of relatively high order in this expansion. 

To simplify the notation, we introduce the di¬ 
mensionless parameter 


c = 


a 

M 2 ■ 


(31) 


As a result of our approximation scheme, we ex¬ 
pand all quantities, such as the metric functions 


^ Note that the constraint = glIM'^ < 0.691 is valid 
for nonspinning solutions at finite coupling. The precise 
value of the upper bound can be modified for large rota¬ 
tion rates m- Indeed, as discussed later in this section, 
the BH mass acquires O(x^) corrections which can be re¬ 
absorbed in the definition of the mass. Nonetheless, the 
bound on afM'^ emerges only from the nonperturbative 
BH solutions and does not appear in the small-coupling 
approximation (to any order in <^). 
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and the scalar field, in terms of the two paramet¬ 
ers C I X- For example, 

= (32) 

Z = 1 j = l 

where the double superscript denotes the or¬ 
der of the expansion in the BH spin parameter 
and in the EDGE coupling parameter, respect¬ 
ively; is the Schwarzschild metric. In prac¬ 

tice, using the spin decomposition previously dis¬ 
cussed, we simply expand the set of radial variables 
/ = {F, as 

/ = EC'/^'^ (33) 

j=0 

where are radial functions which do not de¬ 
pend on the coupling parameter C. By replacing 
these expressions into the field equations derived in 
Sec. |II C[ and solving them order by order in we 
obtain the desired expansion for the metric tensor 
and the scalar field. Remarkably, this procedure 
yields an analytical solution. We compute the ex¬ 
plicit solution up to 0(C^,x®), but the procedure 
can be straightforwardly extended to higher order 
both in C and in y. 

Solving the differential equations at each order 
in X and C yields some integration constants, which 
are uniquely fixed by requiring that 

1 . the metric is asymptotically flat, and the 
scalar field vanishes at spatial infinity; 

2 . there exists an event horizon, where perturb¬ 
ations are regular; 

3. the physical mass and angular momentum of 

the BH are given by M and as meas¬ 

ured by an observer at spatial infinity. In 
particular, the bare mass of the 0 {(^) solu¬ 
tion acquires some corrections to each order 
in (^, which are reabsorbed in the physical 
mass M. 

We note that only one of the two integration con¬ 
stants appearing in the solution of the scalar field 
at each order in C is fixed by requiring regularity 
outside the horizon, while the metric is regular for 
each value of the remaining constants. Although 
this is not evident in the Schwarzschild coordin¬ 
ates adopted here, it can be nonetheless checked 
by computing some curvature invariants. However, 
the remaining integration constants can all be re¬ 
absorbed in the definitions of the physical mass 
and angular momentum, so that the final solution 
truncated at a given order depends only on two 
parameters, as in the Kerr case. 


The explicit expressions of the metric tensor 
and the scalar field up to 0 (('^,x^) are quite long 
and are available in a Mathematica notebook 
provided in the Supplemental Material. For com¬ 
pleteness, the explicit Kerr metric to O(x^) in the 
Hartle-Thorne coordinates is given in Appendix [A| 


III. GEOMETRICAL AND GEODESIC 
PROPERTIES OF THE SOLUTION 


Here we study the properties of the analytical 
solution we have derived. To this aim, we compute 
some geometrical and geodesic quantities which 
characterize the spinning EDGE BH solution to 


A. Event horizon, ergosphere, intrinsic 
curvature, and dilaton charge 


The event horizon is given by the largest root 
r = rh of the equation (cf. e.g., [40]) = 

0 , which yields the following power expansion in 
terms of C and x- 


7 

^ = E + c^x^ + d^x^ + e.iX^ + ftX^) , 

i^O 

(34) 

where the coefficients (a^, bi, Ci, di, Ci, fi) are listecj^ 
in Table |TT| of Appendix]^ As in the Kerr case, the 
horizon radius rh does not depend on the angular 
coordinates. Nonetheless, its intrinsic geometry - 
as computed by considering a spatial section dt = 0 
at r = rh - is nonspherical. Indeed, 

ds^_const,r—Th O^dXt , (35) 


and since ggg explicitly depends on 0, the in¬ 
trinsic geometry is nonspherical. For the line ele¬ 
ment (35), the curvature radius is 


_ cot Og'gg ^ g^ ^ g^ 

gee gle die 9ee ’ 


(36) 


where (only in the above formula) a prime denotes 
differentiation with respect to 9, and for our solu¬ 
tion is 


7 

= 'y ( C[h + X^(™z + Tli COS 26) -|- 
i=0 

+ COS 29 + Ui cos 40)] ; 


(37) 


^ For the sake of clarity, the coefficients shown in the Ap¬ 
pendix will be rounded to some numerical factors. The 
exact expressions are available in the Supplemental Ma¬ 
terial Mathematica notebook. 
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this is constant only when C = 0 = y. Hereafter, 
we adopt the same expansion of Eqs. (34) and (371 
for other physical quantities. The numerical values 
of the coefficients of these expansions are given in 
Appendix]^ whereas their exact form is provided 
in the Supplemental Material Mathematica note¬ 
book. 

The location of the ergosphere is given by the 
largest root of gtt =0, 


the dilaton field ^ ^T) jr and reads 

2 ? J ' 

+ ^ (39) 

where the coefficients di and fi identically vanish. 


B. Moment of inertia 


7 

= '^C[k + +n^cos2e) + 

+ X^{Pi + Qi cos 29 + Ui cos 4:9)], (38) 

where the only nonvanishing spin corrections cor¬ 
respond to even powers of %. 

Finally, the dilaton charge V can be extracted 
from the leading-order large-distance behavior of 


The moment of inertia is defined as / = J/Hh, 
where J is the BH angular momentum and Hh is 
the angular velocity at the horizon of locally non¬ 
rotating observers, 

Hh = — lim . (40) 

In our case we obtain 


^ = 4 - 0.2625000C^ - 0.1721966C^ - 0.1458764C'‘ - 0.1409996C® - 0.1474998C® - 0.1627298C^ 

- _ 0.2359276C^ - 0.2175544C^ - 0.2431079C‘‘ - 0.2776072C® - 0.3283860C® -t 

- 0.3984877C^] -b x‘‘[0-25 - 0.1170266C^ - 0.04956483C^ -b 0.01732049C^ -b 0.09842336C'^ -b 

-b 0.2055222C® -b 0.3503737^] , (41) 


where again the only nonvanishing spin corrections 
correspond to even powers of x- 


C. Quadrupole moment 

According to the BH no-hair theorems, the 
quadrupole moment (as well as the higher-order 
multipole moments) of any regular, stationary, 
asymptotically flat BH in GR is uniquely determ¬ 
ined by its mass M and angular momentum J mi¬ 
lls]- A deformed Kerr geometry as the one just 
discussed, does not necessarily possess this unique 
no-hair property. Since the dilaton charge of this 
solution is not an independent parameter, the mul¬ 
tipole moments of an EDGB BH can all be written 
in terms of M and J, but the relations among them 
will change with respect to Kerr. The ( corrections 
to the BH quadrupole moment are, thus, relevant 
to test the Kerr hypothesis US]. 

To compute the quadrupole moment, we follow 
the general approach described in [44] , in which the 
multipole moments of an asymptotically flat geo¬ 
metry are read off the asymptotic behavior of the 
metric. This approach requires the metric to be 
expressed in asymptotically Gartesian and mass- 


centered (AGMG) coordinates. In particular, in 
order to extract the quadrupole moment, the met¬ 
ric has to be AGMG-2, i.e., gu and {i,j ^ t) 
should not contain any angular dependence up to 
0(l/r^) terms. In our case, the coordinate trans¬ 
formation that enforces such property is 


X 


1 + "- 
r 


2 r 
X^M^ 

9 -\ -— sin 9 cos 


2M2 

j,2 


M(6M - r) 



and does not involve the EDGB coupling C or spin 
corrections higher than second order. In the new 
AGMG-2 coordinates, the gtt component reads 


„ _ 1 ^ rn w20 

gtt--l + — + ^[Q2oy 


+ {I = 0 pole)] -b O 



(42) 


where I 20 is the {I = 2,m = 0) spherical har¬ 
monic, and Q 20 is the m = 0 mass quadrupole 
moment. From our explicit solution, we obtain to 
order 0 {C,X^), 













Q 2 o = - [(1 + 0.1061619C^ + 0.07524246C^ + 0.07459416C'‘ + 0.07756926C® + 

+0.08553316C® + 0.09805643C’^) - (0.0308519 + 0.0408857C + 0.0638894C^ + 

+0.0866408C^ + 0.116314C‘‘ + 0.154763C®)] . (43) 


Interestingly, the 0{x^) corrections to the quadru- 
pole moment are proportional to i-e., they van¬ 
ish in the GR limit. For C ~ 0.4, the 0(C^, X^) cor¬ 
rection to the quadrupole moment relative to the 
Kerr case is about 1.7%, whereas the cor¬ 

rection is approximately 0.5%. Finally, for ( ~ 0.4 
and X ~ 0 . 6 , the 0 ((^,x‘^) correction is approxim¬ 
ately 0 . 1 %. 

We remark that the quadrupole moment of spin¬ 
ning EDGE BHs has been computed numerically 
in [3D]. Our solution has the advantage of giving 
this quantity in analytical form. 


D. Geodesics and epicyclic frequencies 


We shall now consider timelike geodesics in the 
slowly rotating EDGE BH spacetime. We assume 
a minimally coupled test particle and restrict to 
equatorial orbits, for which d = 7r/2 and dO — 0. 
We first compute stable circular orbits; then, by 
considering small perturbations of these orbits, we 
derive the epicyclic frequencies Wr and wg (see, 
e.g.. Refs. [3TH33] b For a stationary-axisymmetric 
spacetime, the ISCO corresponds to the radius at 
which the second derivative of the effective poten¬ 
tial 

V(r) = — ( + 2gLgty + L'^gtt _ 

9rr \ ~ 9ttgipip J 


vanishes. Here, we have introduced the particle- 
specific energy and angular momentum £ and 
L |45|, given by 


£ = - 


L = 




9t(p 4 ” 9ipip^(p 


\J Iht ‘^gtip^ip 9ipip^'^ 


(45) 

(46) 


where Lo^p is the azimuthal angular velocity 


— 


~9tip,r + 





(47) 


Solving V"{r) = 0 order by order, we obtain the 
ISGO radius up to 0(C^,x^): 

+ ^^X + c^X^ + diX^ + 

+ ea^ + fa^). (48) 


Orbits with radius r > risco are stable. Under a 
small perturbation, a massive particle orbiting in 
one of these stable, circular orbits oscillates with 
radial and vertical frequencies given by [3TI - I33] 


2 ^ {gtt +u}p,gt 4 ,y d'^U 

^ 2 grr dr^ 


2 _{gtt+(Apgt^d'^U 

^ ^gee 96“^ 


(49) 

(50) 


These are the epicyclic frequencies. Here U = 
g** — 2 lg*'f’ + with I — LjE being the ra¬ 

tio between the particle angular momentum and 
its energy [^. The full expressions for ujr,ujei as 
well as for oj^ as functions of (r, M, x, C) and up to 
order 0(C^, X^) are available in the Mathematica 
notebook provided as Supplemental Material. We 
explicitly show here their values at the ISGO: 


I ISCO ='^C{ai + hx + CiX^ + d^x^ + 

i^O 

+ SiX^ + fiX^) , (51) 

7 

^‘^^lisco ='^C{ai + hx + CiX^ + d^x^ + 

i=0 

+ f^iX^ + ftX^) , (52) 


whereas uJr \ jg^o = 0 as in the Kerr case. 


E. Comparison with previous results 

As a check, we can compare our results with 
those derived in 122], where the metric of the 
EDGE spinning BH was found to C>(C^,X^) ia 
Boyer-Lindquist coordinates. A direct compar¬ 
ison of the metric coefficients is not possible, since 
the BH solutions have been derived on different 
charts. However, we can overcome this problem 
by computing the Kretschmann invariant JC = 


9<p<p,r 
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Ra 0 'ysR°‘^^^ and evaluating it at a specific point. 


From our solution truncated at we get 


M2 144M2 
JC{r,9)= 48— + ^^ 


M 

(1 — 8cos^ 6 ) 4 -sin^ 9 + 2—^(3cos^ 0—1) 


M® 128 MS 

+14— + — — - 1680 


5 

348909 M^ 
350 7^ ~ 

80334 M^ 
875 7^ 
2168052 Mio 


rp I pC 

304938 M® 


r 

• 2,2 


c' 


2M3 




'M3 

^3 


175 rS 
19638 M® 
' 175 7^ 
59544 Mil 


28023 M^ 
79“ 


35 


r-lO 


35 + 

234816 M® 
1225 7^ 

Mi2 

- 65952 


54431 M^ 

' 1750 7^ 

359468+/!° 

35 Tio” ' 
1448877 Ml' 
rT 


12846 M° 


175 r 
53848 Mil 


M^ 

~^i3 ^ J.4 

77047 M° 
' 1225 7^ 
Mi 2 


144—^ 


„ii 


cos 


350 
■2 0 


5 r 
711114 MS 
' 175 7^ 


- 21984 


„12 


+ 


92679 M° 
35 7^ 


(53) 


Replacing the explicit expression for in Eq. (34) 
we find that on the horizon 


/C(rh,7r/2) = 


+ 


4M4 

C" 


M* 


_ + 

8M4 

327 404023x^ 

127) 784000 


(54) 


This result coincides "with the Kretschmann scalar 
derived in m and evaluated at the event hori¬ 
zon on the equatorial plane in Boyer-Lindquist co¬ 
ordinates. Finally, we have verified the agreement 
between the expression for Moj^ at the ISCO - 
■which is also a gauge invariant quantity - obtained 
from the metric derived in Ref. [27j . and the same 
expr essi on obtained truncating the expression in 
Eq. ® toO{e,X^)- 


F. Accuracy of the expansion 


In this section, we estimate the accuracy of our 
perturbative scheme. In particular, we estimate 
the truncation error arising from neglecting 0 {(^) 
terms in the expansion. To this aim, we compare 
our results with those obtained in Refs. UniUHlIM], 
where a solution for slowly rotating BHs in the 
EDGB theory has been derived at first order in 
X and is “exact” in C (i.e., with no perturbative 
expansion in C). To be con sist en t, w e neg lect t erm s 
ofthe ordev 0{x^) in Eqs. (34), (48), (51) and (52). 


In Fig. [R we compare the dilatonic charge com¬ 
puted in [20l 128] nonperturbatively in <C, with the 
expression in Eq. (39), for y = 0, truncated at 


various orders of C- As expected, for ^ 0.2, 

higher-order corrections are negligible, but they 
contribute significantly as C —t 0.691. To C>(C^)) 
the deviation from the exact result is about 1% 
for ~ 0.6 and is as large as 5% for C 0.691. 


In contrast, the 0{C,'^) truncation differs by about 
30% as c) increases to its maximum value. 



Figure 1. Top panel: Dilatonic charge as a function 
of ^ = a/M^ computed in |20l 128] (g ray markers) 
compared with the expression in Eq. ( 391 truncated at 
0(7), 0(7), and 0(7) and for vanishing spin. Bot¬ 
tom panel: Relative discrepancy between the perturb¬ 
ative and nonperturbative estimates of the dilatonic 
charge, as a function of a for various truncations. 


Likewise, for the set of quantities / = 

{?'h, nsco, Mw^^lisco, Mwelisco}, we compute the 
relative error 


En = - 1 , (55) 

where / represents the exact quantity (nonper¬ 
turbative in 7 [MllMj- We estimate e„ at various 
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orders of approximation in C, for different values 
of the BH spin parameter. In Table |lj we show the 
largest relative errors obtained for all considered 
quantities, at different levels of accuracy, in the 
limiting case ( = 0.691 (left) and ( = 0.576 (right). 
We remark that ( = 0.691 is an extreme situation, 
since for slightly smaller values of ( (i.e., ( = 0.576) 
the deviations are much smaller. 

Figure and Table [l] show that our ana¬ 
lytical solution approaches the exact solution 
of Uni [2H] as the value of n increases, i.e., 
when we consider more and more terms in the 
small-coupling expansion. In particular, for 
nsco,Mu}^\iSco,Mu}e\isco, the relative errors 
(for n = 7) are always smaller than 1% for any 
value of even for the maximum allowed value, 
C 0.691. For the horizon, the threshold above 
which e „=7 > 0.01 is lower, namely () ~ 0.55. 


G. Are spin corrections important? 

The analysis presented in the Sec. mu shows 
that the metric expanded in powers of C, which 
we derived in a closed, analytic form, is a very 
good approximation of the exact numerical result: 
it reproduces the most relevant geodesic quantities 
within 1% for the maximum value C ~ 0.691 and 
within 0.3% for ( ~ 0.576. It is, therefore, justified 
to adopt such higher-order perturbative expansion 
as a starting point to devise strong-held tests of 
gravity. 

In Ref. [34] , we studied the deviations of the azi¬ 
muthal and epicyclic frequencies in a slowly rotat¬ 
ing EDGE BH to hrst order in the spin. However, 
deviations from the Kerr case should increase with 
higher values of the spin. Indeed, as the spin in¬ 
creases, the ISCO gets closer to the horizon, and, 
therefore, observables from orbits near the ISCO 
probe a region of higher curvature, where the de¬ 
viations should be larger. 

In Fig. [^ we conhrm this claim by showing the 
deviations of the horizon and ISCO locations and, 
most important, of the azimuthal frequency uj^p 
and angular epicyclic frequency wg at the ISCO, 
relative to their values computed using the Kerr 
metric approximated at ©(x®), and as functions of 
C and X- For a fixed value of C,, the percentual er¬ 
rors are systematically larger as the spin increases, 
reaching up to 7% for y = 0.6. This large value 
of the spin parameter should be considered as an 
extrapolation. Indeed, our results neglect terms 
of the order ©(y®), which introduce corrections of 
roughly 5% for y ^ 0.6. 

Our perturbative solution is also useful to es¬ 
timate the convergence properties of the expan¬ 
sion. From the coefficients listed in Table [IT] we 
can compute the ratio of the 0(y") and ©(y"”^) 


corrections for a given quantity. For the angular 
epicyclic frequency MwgliscO: this ratio is roughly 
(0.41,0.39,0.37) for n = (3,4,5), in the extreme 
case y ^ 0.5 and C, ^ 0.5. Therefore, the fifth-order 
spin correction is about 20% of the quadratic one. 
Other quantities show a similar behavior. Clearly, 
the convergence improves for smaller values of y, 
whereas it is almost insensitive to the values of the 
EDGE coupling Q. 

Einally, we note that the percentual error of the 
horizon location is almost insensitive to the spin, 
whereas the epicyclic frequencies are much more 
sensitive to this parameter. 


IV. CONCLUDING REMARKS 

With the advent of precision measurements of 
the spectrum of accreting compact objects, it is of 
utmost importance to devise tests of gravity that 
use these measurements to probe the geometry 
near compact objects. To this aim, we have con¬ 
sidered a specific modified theory - namely, EDGE 
gravity - as a case study. This theory has some ap¬ 
pealing theoretical features; for example, it is free 
from instabilities and circumvents the BH no-hair 
theorems. Furthermore, it modifies GR precisely 
in the strong-curvature regions, while passing all 
current solar-system and binary-pulsar tests |1|. 

Spinning BHs in this theory have been studied in 
the past, both numerically UnHUHSO] and analyt¬ 
ically [24U27] to leading order in the coupling para¬ 
meter. Numerical solutions have the advantage of 
being general, but they are impractical for some 
applications, for example for Monte Carlo sim¬ 
ulations spanning a high-dimensional parameter 
space. Approximate analytical solutions can be 
very useful for this purpose, although they are usu¬ 
ally perturbative. 

Here, as a first step to develop precision tests of 
gravity based on geodesic motion near stationary 
BHs, we have constructed an analytical, approxim¬ 
ate solution of EDGE theory describing a deformed 
Kerr BH. The solution is valid to fifth order in the 
spin and to seventh order in the coupling para¬ 
meter, thus, extending previous solutions that are 
valid only to quadratic order in the coupling and 
the spin. With the analytical solution at hand, it 
is straightforward to compute various quantities of 
interest. We have presented the corrections to the 
horizon and ergoregion location, moment of inertia 
and quadrupole moment relative to the Kerr met¬ 
ric, as well as the charge of the dilaton field that 
characterizes this solution. For a given value of the 
coupling C, the solution depends only on the mass 
M and on the dimensionless angular momentum y, 
while the dilaton charge is fixed in terms of M. In 
addition, we have computed some geodesic quant- 
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X 

e„=2(%) 

en=4(%) 

en=6(%) 

en=7(%) 


X 

£„=2(%) 

e„=4(%) 

en=6(%) 

en = 7(%) 

rh/M 

0 

5.90 

4.45 

3.72 

3.48 

rh/M 

0 

1.33 

0.52 

0.24 

0.17 

nsco/M 

0 

1.00 

0.58 

0.43 

0.39 

nsco/M 

0 

0.32 

0.093 

0.038 

0.026 


0.05 

1.11 

0.65 

0.49 

0.44 


0.05 

0.37 

0.12 

0.055 

0.042 


0.10 

1.23 

0.72 

0.54 

0.49 


0.1 

0.42 

0.14 

0.074 

0.059 

Mu}^p\isco 

0 

1.36 

0.79 

0.59 

0.53 

Muj^p\isco 

0 

0.44 

0.13 

0.053 

0.036 


0.05 

1.56 

0.95 

0.72 

0.66 


0.05 

0.56 

0.23 

0.14 

0.13 


0.10 

1.88 

1.22 

0.98 

0.91 


0.1 

0.80 

0.44 

0.35 

0.33 

Muie\isco 

0.05 

1.53 

0.92 

0.69 

0.63 

Mlub\isco 

0.05 

0.54 

0.20 

0.12 

0.10 


0.10 

1.78 

1.13 

0.88 

0.81 


0.1 

0.71 

0.36 

0.27 

0.25 


Table 1. Left: the relative error e„ [cf. Eq. (55l] between different quantities listed in the first column, computed 
through the solution derived in |28| . nonperturbative in and compared with our perturbative results truncated 
at 0{C,'^). We consider the maximum value of C, allowed for BH solutions in EDGE gravity, = 0.691, and 
different values of the BH spin parameter. Right: Same for C, = 0.576. 




a//W^ a//W^ 


Figure 2. L eft p anel: Th e percentual error in the horizon radius rh and in the ISCO nsco for our perturbative 
result [Eqs. ( |34[ | and (481], relative to the Kerr solution expanded at 0{x^), as a function of the EDGE coupling 
parameter a. We consider three values of the BH spin parameter x = (0.2,0.4,0.6). Right panel: Same as the 
left panel but for the epicyclic frequencies Eqs. (|51ll and (1521) evaluated at the ISCO. 


ities, namely, the ISCO location and the azimuthal 
and epicyclic frequencies as functions of M, x, 
and the orbital radius r. 

When truncated at first order in the BH spin, 
our solution reproduces the most relevant geodesic 
quantities obtained in [2Q1 [28] with a numerical 
approach, within 1% for the maximum value ^ 
0.691, and within 0.3% for C, ~ 0.576. The accur¬ 
acy of the solution grows dramatically for smaller 
values of the coupling. These results indicate that 
our perturbative solution is a good approximation 
of the exact numerical results. 

In a future publication, we will extend the ana¬ 
lysis of Ref. |Mj, which studied how observations 
of quasiperiodic oscillations in the spectrum of ac¬ 


creting BHs can be used to constrain EDGB the¬ 
ory in the strong-held regime, to larger values of 
the BH spin. A similar analysis can also be per¬ 
formed for other tests based on stationary BHs, for 
example, those based on the broadened iron line 
(e.g., [Ml [47]) or the continuum htting method 
(e.g., [Ml; see also, [35|)- On the technical side, 
our perturbative approach is generic: it can be ap¬ 
plied to any order in C and in x, as well as to 
other modihed-gravity theories and different spin¬ 
ning solutions. 
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Appendix A: Kerr metric in the 
Hartle-Thorne approximation 

In this appendix, we show the form of the BH 
solution in the Hartle-Thorne approximation for 
a = 0, i.e., the slowly rotating Kerr BH in GR, 
up to the fifth order in the BH spin angular mo¬ 
mentum, 


9tt = 


= 1 - 


2M 


6 

5r® 


+ r 


146 


4M\ 


Mr3 

16 


44M 


1 P 2 (cos 0 ) + 


2 cos 29 
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12M 


11 


7r8 7M2r6 7r^ IMr"^ 


24 

+ FN ) cos(2d) -|-sin^(0) 


46 74 
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Appendix B: Coefficients of the small coupling 

In the following table, we show the numerical 
coefficients of various analytic expansions presen¬ 


ted in the previous sections, as functions of the BH 
spin angular momentum and the EDGB coupling 
parameter. For the sake of clarity all the coeffi¬ 
cients are rounded to the seventh digit. The exact 
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expressions are available in a Supplemental Mater- ial Mathematica notebook. 
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Table II. Numerical values of the coefficients of the expressions (341,(391,(481, (51 l-(52l (left panel), and of the 
ergosphere and intrinsic curvature radius (38l], (371 (right panel). 
























